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Abstract 



Dynamical control of biological systems is often restricted by the practical con- 
straint of unidirectional parameter perturbations. We show that such a restriction 
introduces surprising complexity to the stability of one-dimensional map systems 
and can actually improve controllability. We present experimental cardiac control 
results that support these analyses. Finally, we develop new control algorithms that 
exploit the structure of the restricted-control stability zones to automatically adapt 
the control feedback parameter and thereby achieve improved robustness to noise 
and drifting system parameters. 
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L INTRODUCTION 

Recent success controlling complex dynamics of nonlinear physical and chemical sys- 
tems [p]-|T^|] has opened the door for the control of biological rhythms. Some researchers 
have speculated about the medical implications of controlling heart-rate dynamics or brain 
rhythms [|T7]-^2p. However, biological systems typically have characteristics that require special 
consideration. For example, biological control studies to date [ lT7|^rrS|P^| - P5|PT| , p7| -PD|] have 
required that the control interventions be unidirectional — only allowing shortening of a param- 
eter. Such a restriction is somewhat analogous to trying to balance a broomstick vertically on 
one's palm using horizontal hand movements in only one direction. Intuitively, one might expect 
that such a restriction would limit controllability. However, as we will demonstrate in this paper, 
such a restriction introduces some surprising complexity to the stability properties of controlled 
one-dimensional map systems. 

In fact, the unidirectional restriction can actually improve the controllability of some sys- 
tems [BT]]. In this paper, we will show how restricted control can introduce stability zones that do 



not exist in the unrestricted case. Furthermore, we will show that some of these zones were present 
in recent cardiac control experiments [^TJ. Finally, we will exploit the structure of the stability 
zones in robust new control algorithms that automatically adapt the control feedback parameter. 



H. DELAYED FEEDBACK CONTROL OF SYSTEMS DESCRIBED BY ONE-DIMENSIONAL 

MAPS 

In this study we will consider the control of systems whose dynamics can be described by 
one-dimensional maps: 

X n+1 = f(X n ,\), (1) 

where X n is the variable to be controlled and A is an experimentally accessible system parameter. 
The goal is to stabilize the system state point = [X n , X n -i] about an unstable period- 1 fixed 
point £* = [X*, X*], where X* = f(X*, A), by perturbing A by an amount: 
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5\ n = —{X n -i — X n ), (2) 

where a is the feedback gain parameter. 

The advantage of such a control scheme is that relatively little a priori system information is 
required to implement control and stabilize £*. The only requirement is knowledge of the sign of 
|j so that perturbations can be applied in the correct direction. In fact, knowledge of the value of 
£* is unnecessary because the controlled system's fixed point is identical to that of the uncontrolled 
system. Furthermore, if the fixed point drifts during the course of the control (as is common for 
biological systems), the controlled system will track the fixed point, provided that the system stays 
in the stable range of the feedback gain parameter a. 

The control algorithm of Eq. ^ is an example of delayed feedback control, a technique that 
has been used in a variety of modeling and experimental studies [|2|-p], |32|JlO| , ^4| ]. In section |V 



we 



will present an example of a biological system with constraints that restrict the control algorithm 
— only allowing unidirectional perturbations of A. The purpose of this study is to examine the 
implications of such a restriction. 

A. Linear stability analysis of unrestricted delayed feedback control 

For unrestricted control, in which 5X n can be positive or negative, linearizing the controlled 
system about a fixed point at the origin gives: 

X n+1 = AX n + P(Y n -X n ), (3) 
Y n +i = X n , 

where A=fe\e and /? = f(§)| e - 



The eigenvalues of Eq. [| are (A — (3 ± (A — j3) 2 + 4(3) /2. The fixed point is stable provided 
that both eigenvalues fall inside the unit circle. This condition is met when: 

-1<P<~(A + 1), (4) 



Hall and Christini 4 

for A < 1 J21|j33| ] . Note that the stability zone shrinks to zero for A < —3, thereby limiting 
the applicability of the unrestricted control algorithm to maps with a sufficiently shallow slope 
(i.e., — 3 < A < 1) at £*. Furthermore, for A > 1 there exists no real value for j3 such that 
the eigenvalues fall within the unit circle. Thus, unstable positively-sloped fixed points cannot be 
stabilized. 



III. RESTRICTED DELAYED FEEDBACK CONTROL 

Restricting the above control algorithm by only allowing shortening of A gives the following 
controller: 

5\ n = e n ^(X n „ 1 -X n ), (5) 

where 

(i if(x n -r n )>o, 
e n = I (6) 

otherwise. 

Thus, when n = 1 the control is active (i.e., a perturbation is delivered), and when n = the 
control is inactive (i.e., no perturbation is given). 

A. Linear stability analysis of restricted delayed feedback control 

The restricted control algorithm of Eq. |5] gives the following linearized controlled system: 

X n+1 = AX n + Q n (3(Y n - X n ), (7) 

Geometrically, the restriction of Eq. ^| means that perturbations will only be applied if the state 
point lies above the return-map line of identity X n+ i = X n . The dynamical effects of this 
restriction depend on the sign of the slope at £*. 
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1. Restricted control for A < — 1 

Typical dynamics of restricted control for negatively-sloped unstable directions (A < — 1) are 
depicted in Fig [I]. This figure shows eight control trials of a linear map with A = — 4 for different 
values of [3; there are four examples of stable control and four examples of unstable control. 

Figure []]( a) shows a case in which the restricted control algorithm failed to stabilize the unstable 
fixed point £* [which is located at the intersection of the uncontrolled system map (solid line) and 
the line of identity (dotted line), and is denoted by a solid triangle] with f3 = —2.80. The dot-dash 
lines correspond to the system maps when A is perturbed. A series of arrows originate at the initial 
state point and connect consecutive state points (solid circles, numbered consecutively). In this 
case, the initial state point (1) is followed by a control intervention, which causes the next iterate 
(2) to fall below the line of identity. According to Eq. the next iterate (3) will be uncontrolled 
and therefore will fall on the solid line. Furthermore, because the first controlled iterate (2) was 
less than the fixed point, the next iterate (3) will be above the line of identity, leading to a control 
intervention at the fourth iterate (4). Thus, control is applied every other iterate so that the sequence 
of 0„ is 0101.... In this case, the fixed point is not stabilized because control fails to direct the 
system closer to the fixed point (i.e., the arrows spiral away from £*). 

Figure [T](b) shows control with the same value of A, but using a slightly more negative [3 
value, f3 = —3.1. For these parameter settings, it can be seen that control is also applied every 
other iterate so that the sequence of 6„ is again 0101.... However in this case, the fixed point 
is stabilized successfully because the control interventions direct the system closer to the fixed 
point (i.e., the arrows spiral towards £*). Note that in order to maximize the clarity of the control 
sequence diagram, the axes in panel (a) are not scaled the same as those for panel (b). Similarly, 
for all of Fig. [IJ, axes from different panels are not necessarily scaled the same. 

Figure |T](c) shows a trial in which [3 was decreased to —3.23. As in the previous example, 
the first controlled iterate (3) is below the line of identity (dictating that the next iterate (4) is un- 
controlled). However, in this case the perturbation is larger than would occur with the parameter 
settings of Fig. [IJb), such that the controlled iterate (3) is slightly larger than the fixed point. This 
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dictates that the next iterate (4) is below the line of identity, which leads to a second consecutive 
uncontrolled iterate (5). Thus, control is applied in a 001001... sequence. This sequence is sta- 
ble for (3 = —3.23 because the control perturbations direct the system closer to the fixed point. 
However, if (3 is decreased to —3.40, it can be seen from Fig. [I](d) that this generates a 001001... 
control sequence that is unstable because the system is directed away from the fixed point. 

If (3 is made more negative, then a new control sequence is achieved. Figure [l|(e) shows an 
unstable 011011... sequence for (3 = —5.50. In this case, the first perturbation is so large that the 
controlled iterate (2) is above the line of identity, dictating that the next iterate (3) is also controlled. 
The second controlled iterate (3) is below the line of identity and below the fixed point, thereby 
producing the 01 101 1 ... sequence. In this case, the fixed point is not stabilized because the control 
perturbations moved the state point away from the fixed point. However, when (3 is decreased to 
—5.76, as seen in Fig. |T](f), a 01 1011... sequence stabilizes the fixed point. 

Like the transition from the stable 0101... sequence to the stable 001001... sequence depicted in 
Figs. [j](b) and (c), there is a transition from the stable 01 101 1 ... sequence to the stable 001 1001 1... 
sequence as (3 is decreased further. Figure |T](g) shows the stable 001 1001 1... case for (3 = —5.798. 
The 00110011... sequence becomes unstable as (3 is decreased still further; Figure [T|(h) depicts 
this case for (3 = —5.95. 

The progression of unstable and stable periodic control sequences continues indefinitely as (3 
is decreased. In fact, the switching parameter 0„ imposes the following progression of control se- 
quences as (3 is decreased from zero: unstable 01 1 , stable 01 1 , stable 001 1 , unstable 001 1 , unstable 
01 2 , stable 01 2 , stable 001 2 , unstable 001 2 , unstable 01°°, stable 01°°, stable 001°°, unstable 
00 1°°, where l fc denotes k consecutive ones (control perturbations) before the sequence repeats 
itself. 

Because of the progression of the control sequences imposed by the switching term n , X k+ i 
can be expressed as: 

Xk+i = efcX , (8) 
where e& is given by the following iterative expression: 
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e k = (A- /3)e fc _! + (3e k _ 2l (9) 

with eo = A for all sequences and e± = A 2 + (3(1 — A) for the 01 fc sequences or e\ = A 2 for the 
001 k sequences. 

The boundaries of the stability zones can be computed by using the criterion that stable se- 
quences move the system closer to the fixed point after one control sequence. Because X k+1 is the 
last iterate of the first 01 fc control sequence and X k+2 is the last iterate of the first 001 fc sequence, 
the stability conditions are e k < 1 and e k+ i < 1 for the 01 fc and 001 fe sequences, respectively. 
Therefore, the boundaries are given by k degree polynomials in (3. For example, the k = 1 control 
sequences are stable for 1 + A + 1 / A < (3 < 1 + A for A < — 1 [0 . Figure |2| depicts the stability 
zones (shaded regions) for k = 1 and k = 2. The boundaries between the stable 01 fc and 001 fc 
sequences are defined by the condition e k = (dotted curves in Fig. §). These curves mark the 
optimal parameter values for a given stability region, because the fixed point is reached after a 
single control sequence 01 fc . 

The striking feature of this analysis is that the domain of control is extended by the restriction 
of Eq. H [pT]|. In fact, Fig. || shows that for all A < —1, stable control sequences exist for 
the restricted system. This is in contrast with the unrestricted system, for which stable control 
sequences exist only for A > —3, as shown by the dashed triangular region marking its stability 
zone (according to Eq. g). 

While there are an infinite number of stable zones corresponding to an arbitrary number k con- 
secutive control perturbations, the stability zones are bounded by the curve (3 = A — 2 — 2y/l — A. 
This boundary is computed by recognizing that as k approaches infinity, control is always active 
and X n+ i > X n for every iterate. Thus, the algorithm behaves just like the unrestricted control of 
Eq. |] with real eigenvalues greater than one — a condition met only when (3 is below the boundary. 

2. Restricted control for A > 1 

Typical dynamics of restricted control for positively-sloped unstable directions (A > 1) are 
depicted in Fig ^ This figure shows four control trials of a linear map, with A = 2.1, for different 
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values of (3; there are three examples of unstable control and one example where £* is controlled. 

Figure |](a) shows an unstable 01°° sequence for (3 = 1.5. In this case, the perturbations are so 
small that all state points lie above the line of identity. The control slows the exponential growth 
(which would be marked by a rapid exit from £* along the solid line), but fails to force an approach 
to£*. 

If (3 is decreased to — 2.5, then the perturbations are large enough so that the first controlled 
X n is smaller than the previous X„_i. £ n then falls below the line of identity, thereby generating 
the 01 1 control sequence depicted in Fig. |3](b). In this case, the sequence is unstable because a 
given controlled point is further from than the previous controlled point. 

For (3 = 3.5, control is successful; the state point approaches £* in a 01 1 sequence as shown 
in Fig. |](c). Such control is successful for a noise-free model system. However, for a real-world 
system, once the state point is sufficiently close to £*, noise will eventually kick £ n to the opposite 
side of £*. Subsequent £ n will fall below the line of identity, causing the control to be deactivated 
and leading to an exponential departure from £*. 

When (3 is increased further, the first perturbation can be so large that the state point will be 
kicked to the other side of £* as shown in Fig. ^|(d) for (3 = 4.5. Again, control is subsequently 
deactivated and the system diverges from the fixed point. 

The boundaries between the different control sequences are depicted in Fig. |} The unstable 
01°° sequence occurs for (3 < A. This boundary was computed by finding the value for (3 such that 
the first controlled iterate falls on the line of identity; this occurs for (3 = A. For (3 > A, the first 
controlled iterate lies below the line of identity, thereby shutting the control off for the next iterate. 
Thus, (3 = A marks the boundary between the 01°° and 01 1 sequences. The boundary between 
unstable and stable 01 1 sequences was computed by finding the value for (3 such that the iterate 
subsequent to control is equal to the previous uncontrolled iterate; this occurs for (3 — 1 + A. 
For slightly larger (3 values, the iterate subsequent to control is closer to £* than the previous 
uncontrolled iterate. Thus, the unstable 01 1 sequence is bounded by A < [3 < 1 + A. The final 
boundary marks the end of the converging 01 1 sequence and can be found by determining the 
value for (3 such that the first controlled iterate lands at £*. This occurs for (3 = A 2 /{A — 1). 
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For (3 > A 2 /(A — 1), the first controlled iterate, and all subsequent iterates, lie below the line 
of identity thereby shutting of the control. Therefore, the converging 01 1 sequence occurs in the 
shaded region 1 + A < [3 < A 2 /{A — 1) pfJ3]], and the unstable 010°° sequence occurs for 



P > A 2 /(A- I). 

Thus, for A > 1 the best that the restricted control algorithm can offer is a temporary reversal 
of divergence from the fixed point. However, in section |V A| we will make a simple modification 
to the restricted control algorithm for A > 1 that will keep the system in the vicinity of £* in the 
presence of noise. 

IV. EXPERIMENTAL OBSERVATION OF RESTRICTED CONTROL SEQUENCES 



As described in Ref. []2I{], we have studied the control of a particular cardiac conduction in- 
terval, known as the atrioventricular (AV) nodal conduction time, in in vitro rabbit heart exper- 
iments. Because of the nonlinear excitation properties of AV-nodal tissue, the dynamics of AV- 
nodal conduction can bifurcate from a period- 1 regime (where every impulse propagates through 
the AV-node at the same rate) to a period-2 regime (where propagation time alternates in a long, 
short, long, etc. pattern on a beat-to-beat basis) during rapid atrial excitation. It has been demon- 



strated [ p5| , [36| ] that these dynamics can be described by a period-doubling bifurcation of a one- 
dimensional map of the form of Eq. [j] where X n is the AV-nodal conduction time and A n is the 
time between when the AV-node finishes conducting one impulse and when it starts conducting 
the next. 



The goal in Ref. Q2T[ ] was to eliminate the alternating rhythm by stabilizing the underlying 
period one fixed point X*. This was achieved by delivering electrical stimuli to the atrial tissue in 
order to transiently shorten A n . Because there is no practical way to lengthen A n , the timing of the 
electrical stimuli was determined by the restricted controller of Eq. |5|. 

Although the in vitro rabbit cardiac system of Ref. [fTl]] was not linear, application of the 
restricted control algorithm resulted in several of the control sequences predicted in the above 
linear system for A < —1. These control sequences were especially clear at the initiation of 
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control when perturbations were largest. 

For example, Fig. |5](a) shows the variable X n and the control parameter A n during an unstable 
01 1 sequence for a feedback gain a = 3.3 (corresponding to a negative f3 because |£ < in the 
cardiac control experiments). The first controlled beat is indicated by the arrow and corresponds 
to a negative perturbation of Ao (all control perturbations are negative as imposed by the switching 
term n ). Because the system was nonlinear, oscillatory growth of X n was quenched and the 
original large amplitude alternation of X n was reduced in magnitude — but not eliminated. 

When a was increased to 5.0 (as shown in Fig. |](b); corresponding to a later segment of the 
same control trial that is shown in Fig. |](a)), the system shifted to a stable 00 1 1 sequence that 
eliminated the alternation of X n . After the fourth perturbation to A n (beat 303), the system shifted 
to a stable 01 1 sequence. This shift resulted from the close proximity of the 01 1 and 00 1 1 stable 
zones (Fig. Q); noise or drift in the system can cause such transitions. 

Figure |5](c) shows a stable 001 2 control sequence that eliminated the alternation of X n in a 
different rabbit heart using a = 2.5. (Note that the a values from distinct trials are independent.) 
Similar to the sequence transitions in Fig. |5|(b), the system switched to its adjacent stable 01 2 
control sequence shortly after the control was initiated, and later switched back to the stable 00 1 2 
control sequence. 

V. MODIFICATIONS TO THE RESTRICTED CONTROL ALGORITHM 

A. Automatic adaptation of the feedback gain for A < - 1 

Figure §(a) (unsuccessful) and (b) (successful) demonstrate that successful control is depen- 
dent on the proper choice of a. Such dependence is a critical limitation given that the information 
required to determine the correct value of a (A, £*, and cannot be easily determined prior to 

control. Furthermore, the nonstationarities typical of biological systems imply that the appropriate 
value of a may drift over time, thereby increasing the likelihood of control failure if a is fixed. 
To eliminate the limitations of a fixed a value chosen prior to control, we have developed a new 
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technique that adaptively estimates a J37|]. This adaptive approach is especially appropriate for 
applications (e.g., cardiac arrhythmia control) that cannot afford control failure of the type shown 
in the control attempt of Fig. |](a). 

This new technique exploits the structure of the stability zones to automatically adapt a such 
that £* is stabilized more robustly. Because multiple perturbations away from the fixed point are 
not desirable, the optimal stability zone is the k = 1 zone. Furthermore, because the stable k = 1 
zone has the largest area, it will be the most robust to noise and drifting system parameters. To 
target this zone, a is adapted according to: 

+ 5a if e n _ 4 ... 0„_i = 0101 or 1010, 

(10) 

a n -i — 5a otherwise, 

where 5a is a small increment. When (ff)|^* is negative (as in the cardiac experiments of 
Ref. pl|]), a and 5a are positive. Otherwise they are negative. 

The algorithm of Eq. [R] is motivated by examining the stability zones in Fig. ^. For k — 1, 
optimal control occurs when (3 is at the boundary between stable 01 1 and stable 001 1 (dotted curve 
of Fig. 0). If (3 is too large, the control sequence will be 01 1 (unstable if (3 is so large that it is 
above the k=l stability region or stable if (3 is within the stability region but above the optimal 
control boundary). In such a case, the adaptation of Eq. [T0| will decrease (3. In contrast, if (3 is too 
small, the control sequence will be 001 1 or some sequence with k > 1 (unstable if [3 is so small 
that it is below the k = 1 stability region or stable if (3 is within the stability region but below 
the optimal control boundary). In such a case, the adaptation of Eq. [H] will increase [3. Thus, 
the adaptation will adjust the system so that it oscillates between the stable 01 1 and stable 001 1 
sequences, provided that the increment 5a is small enough so that the stepsize of (3 is sufficiently 
less than the height of the k = 1 stability zone. Specifically, the condition \5a\ < 
ensures that the stepsize is less than half the height of the zone. 

To illustrate the adaptive algorithm, we implemented the restricted controller of Eq. ^] with the 
feedback gain a replaced by a n given by Eq. [HJ. a was randomly chosen between —5 and —10 
and 5a = —0.1. We applied this controller to the quadratic map: 
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X n+ \ — A n X n (l — X n ) + ( n , (11) 

where ( n is a normally-distributed random variable with a mean of zero and a variance of 0.001. 
The goal was to stabilize the fixed point X* = (A — 1)/A . Note that for the quadratic map, 
(§01?* = TP^ * s positive. Thus, the sign of the perturbations is opposite to that of the cardiac 
experiments of Ref . [|2~T|] . 



Figure || shows the results of an adaptive control trial of the quadratic map. For 1 < n < 500, 
Ao = 3.30, corresponding to an uncontrolled stable period-2 orbit. Control was initiated at iterate 
n = 125 with an arbitrary initial value of a n = —5.25. The adaptive algorithm rapidly stabilized 
£* as the adaptations of a n kept the system in the k = 1 stability zone. The control was deactivated 
at n = 375 and the period-2 cycle returned. 

At n = 500, A was switched to 3.52, which corresponds to a period-4 rhythm. Control 
was reactivated at n = 625 with an arbitrary initial value of a n = —8.85. The adaptive control 
stabilized £* until the control was turned off at iterate 875. At n = 1000, A was switched to 3.65, 
which is in the chaotic regime. Control was reactivated at n = 1125 with an arbitrary initial value 
of a n = —5.63. £* was stabilized after approximately 160 iterates. When control was deactivated 
at n = 1375 the chaos resumed. The oscillations in a n , as governed by Eq. [TIJ are apparent in 
each control phase of the trial. 

To demonstrate the ability of the adaptive algorithm to track a drifting fixed point, we applied 
the control algorithm to a quadratic map with A = 3.0 and A n increased by an increment of 0.001 
every iterate (A n+ i = A n + 0.001). As seen in Fig. [7], the small increments to Ao introduced a slow 
drift in the system dynamics and fixed point. Control was initiated at iterate 250 while the system 
was in its period-2 regime. Control was maintained for 500 iterates. During this period, X* drifted 
from X* = 0.692 to X* = 0.736. Nevertheless, the control algorithm had no trouble tracking X*. 
In fact, it can be seen that when control was deactivated at n = 750, the system had passed into 
the chaotic regime, an occurrence that did not disrupt control. However, if the feedback gain was 
held fixed rather than adapted, then control could not have been maintained for the entire control 
period (not shown). 
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B. Control of fixed points for A > 1 



In section |III A 2| , we demonstrated that the restricted control algorithm can induce a transient 
approach towards £* when A > 1. However, as mentioned, if the system is kicked to the other 
side of £*, control is deactivated and the system diverges from £*. If the algorithm could detect 
such an occurrence and reverse the sign of the perturbations, then could be approached from the 
opposite side of [|3"%|]. This idea motivates the following modification of the restricted control 
algorithm: 



5\ n = B n — (X n _i — X n ) 



(12) 



where 



e, 



i if$ n (x„-y n ) >o, 

otherwise, 



(13) 



and 



i ife n _ 4 ...e n _ 1 = iooo, 



1 otherwise. 



(14) 



For stable control, a fixed value of the feedback gain (a n = a ) is chosen so that 1 + A < (3 < 



A /(A — 1). In order to adaptively control fixed points with A > 1, the controller in Eq. [12| can 



be used with a modified adaptive feedback gain algorithm: 



a, 



a n _i + 5a if 6 n _ 4 ... @ n -i = 0101 or 1010, 
Oi n -i — $ a otherwise. 



(15) 



Such a combination is feasible because the control-sequence boundaries for A > 1 dictate that 
Eq. [15] will direct the system towards the converging 0101... sequence (see Fig. ||). However, as in 
the case when A < — 1, the increment in the feedback gain should be chosen such that a n remains 
in the stable 0101... zone. Choosing \5a\ < \{A — 1)(|^)|^»| _1 ensures that the increment is less 
than half of the height of the zone. 
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In order to illustrate the control of an unstable fixed point with A > 1, we applied the modified 
control algorithm to the cubic map: 

X n+1 = -4(m + 1)X* + 6(m + l)X 2 n - (2m + 3)X n + \ n + ( n , (16) 

where A n is perturbed according to Eq. [12| with A = 1, m is the slope of the map at the fixed point 
X* = 0.5, and ( n is a normally-distributed random variable with a mean of zero and a variance of 
0.001. 

Figure|| illustrates control of the cubic map without adaptation of the feedback gain (a n is fixed 
at a = 8.0). For 1 < n < 500, m = 2.2, corresponding to an uncontrolled stable period-2 orbit. 
After control was initiated at n = 125, X* was stabilized after about 80 iterates. Stabilization of 
X* was maintained until control was deactivated at n = 375, after which the system returned to 
the period-2 orbit. At n = 500, m was set to 2.7, which moved the system into the chaotic regime. 
Control was initiated at n = 625. After about 10 iterates, X* was controlled for about 20 iterates. 
The system subsequently escaped control for about 30 iterates before control was recaptured and 
maintained until the algorithm was deactivated at n = 875. 

The adaptive feedback gain algorithm is illustrated in Fig. |9|, which shows control of a drifting 
cubic map with m n+1 = m n + 0.001 and m = 2.0. Control was activated for 250 < n < 750. 
The fixed point was controlled successfully during this period. The fixed point location does not 
change for the drifting cubic map, but the system drifts into the chaotic regime by the end of the 
control period. 

VI. CONCLUSION 

Surprisingly, the typical biological restriction of unidirectional control perturbations enhances 
the controllability of fixed points with A < — 1 in systems described by one-dimensional maps. 
Because biological systems typically drift over time, dynamic control algorithms must adapt to 
system nonstationarities. Although the restricted delayed feedback control technique allows for 
moderate tracking of the fixed point as long as the system remains within a stability zone, it is ill- 
suited for systems with significant drift. For such systems, automatic adaptation of the feedback 
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gain parameter ensures that the drifting system is directed to, and remains within, the largest stabil- 
ity zone. Thus, with the dual benefits of the increased stability of unidirectional restricted control 
and the adaptability of on-the-fly gain estimation, such control techniques could be of significant 
value to the control of biological systems. Indeed, a recent set of clinical experiments [^yp have 
shown that adaptive restricted control of this type can successfully eliminate the same alternating 
AV-nodal conduction rhythm that was controlled in the rabbit experiments of Ref . ^TJ] . Further- 
more, we have shown that simple modifications of the restricted control algorithm can control 
fixed points with A > 1 — an impossible task for the unrestricted feedback controller. Thus, this 
algorithm may also have applicability in physical systems that allow bidirectional perturbations. 
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FIGURES 

FIG. 1. Return maps showing the progression of control sequences as A is decreased below —1. Se- 
quential state points are numbered, the dotted diagonal line is the identity line X n+ \ = X n , the solid line is 
the map of the uncontrolled system with slope A = —4, the fixed point £* is denoted by the solid triangle, 
and the dot-dash lines show the system map when perturbed by control interventions, (a) f3 = — 2.8 results 
in an unstable 01 1 control sequence, (b) (3 = -3.1; stable 01 1 . (c) (3 = -3.23; stable 001 1 . (d) (3 = -3.4; 
unstable 001 1 . (e) (3 = -5.5; unstable 01 2 . (f) /3 = -5.76; stable 01 2 . (g) (3 = -5.798; stable 001 2 . (h) 
(3 = —5.95; unstable 001 2 . Note that axes from different panels are not necessarily scaled the same. 

FIG. 2. Stability zones of unrestricted (Eq. |3|) and restricted (Eq. [Sj> control for A < — 1. The triangular 
region enclosed by the dashed lines in the upper right corner is the stability zone for unrestricted control. For 
restricted control, the k = 1 (01 1 and 001 1 ) and k = 2 (01 2 and 001 2 ) stability zones are the shaded regions 
enclosed by the solid curves (which are A;-degree polynomials in (3, determined via Eq. as described in the 
text). The dotted curves inside the zones mark the transition from 01 fe to 001 fc . The annotated open circles 
a-f correspond to the control parameters for panels (a)-(f) of Fig. [I[ The three vertically-spaced solid dots 
indicate that there are an infinite number of stability zones for larger k. The infinite sequence of stability 
zones is bounded by the curve f3 = A — 2 — 2\/I — A. 

FIG. 3. Return maps showing the progression of control sequences as A is increased above 1. Sequen- 
tial state points are numbered, the dotted diagonal line is the identity line X n+ i = X n , the solid line is the 
map of the uncontrolled system with slope A = 2.1, the fixed point is denoted by the solid triangle, and 
the dot-dash lines show the system map when perturbed by control interventions, (a) /3 = 1.5 results in an 
unstable 01°° control sequence, (b) (3 = 2.5; unstable 01 1 sequence, (c) (3 = 3.5; converging 01 1 sequence, 
(d) (3 = 4.5; unstable 010°° sequence. Note that axes from different panels are not necessarily scaled the 
same. 
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FIG. 4. Stability zones of restricted control (Eq. ||) for A > 1. The zone of semi-stability (shaded 
region; denoted 01 *) and the different zones of instability are separated by the curves (3 = A, (3 = 1 + A, 
and j3 = A 2 /(A — 1). The annotated open circles a-d correspond to the control parameters for panels 
(a)-(d) of Fig. | 

FIG. 5. Control sequences observed in the rabbit heart experiments of Ref. [21]. The first control 
perturbation in each panel is indicated by an arrow, (a) An unstable 01 1 sequence for a = 3.3. (b) The 
same preparation with a = 5.0. In this case the control begins in a stable 001 1 sequence and shifts to a 
stable 01 1 sequence, (c) A stable 001 2 sequence, shifting to a stable 01 2 sequence, and returning to a stable 
001 2 sequence in a different preparation with a = 2.5. 



FIG. 6. Adaptive control of the quadratic map of Eq. [11|. Control was activated during the intervals 
labelled with a "C". Iterates 125-375, 625-875, and 1125-1375 show control of period-2 (A = 3.30), 
period-4 (Aq = 3.52), and chaotic dynamics (Ao = 3.65), respectively. 



FIG. 7. Adaptive control of a drifting quadratic map (Eq. [llj). The baseline value of Ao = 3.00 was 
incremented by 0.001 each iterate (A n +i = A n + 0.001), causing a slow drift in the system. Control was 
activated from 250 < n < 750 (labelled with a "C"). During this period, X* drifted from X* = 0.692 to 
X* = 0.736. The adaptive algorithm tracked the drifting fixed point as the system moved into the chaotic 
regime. 



FIG. 8. Control of the cubic map of Eq. |16j with with Ao = 1.0. Control was activated during the 
intervals labelled with a "C". Iterates 125-375 and 625-875 show control of period-2 (m = 2.2) and chaotic 
dynamics (m = 2.7), respectively. 



FIG. 9. Adaptive control of a drifting cubic map (Eq. 16). The baseline value of itlq = 2.0 was 
incremented by 0.001 each iterate (m n+ i = m n + 0.001), causing a slow drift in the system. Control was 
activated from 250 < n < 750 (labelled with a "C")- The fixed point location does not change for the 
drifting cubic map, but the system drifted into the chaotic regime by the end of the control period. 
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